Extraction of optical Bloch modes in a photonic-crystal waveguide 
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We perform phase-sensitive near- field scanning optical microscopy on photonic-crystal waveguides. 
The observed intricate field patterns are analyzed by spatial Fourier transformations, revealing sev- 
eral guided TE- and TM-like modes. Using the reconstruction algorithm proposed by Ha, et al. 
(Opt. Lett. 34 (2009)), we decompose the measured two-dimensional field pattern in a superposi- 
tion of propagating Bloch modes. This opens new possibilities to study specific modes in near-field 
measurements. We apply the method to study the transverse behavior of a guided TE-like mode, 
where the mode extends deeper in the surrounding photonic crystal when the band edge is ap- 
proached. 



INTRODUCTION 

Near-field scanning optical microscopy (NSOM) is a 
powerful tool to study objects with a resolution below 
the diffraction limit pQ. A unique feature of NSOM is 
the ability to tap light from structures that are designed 
to confine light, such as integrated optical waveguides 
[2j [3] and cavities [3H7] . Using the effect of frustrated to- 
tal internal reflection, light that is invisible to other mi- 
croscopy techniques can be detected. It is mainly for this 
reason that NSOM is so useful in the study of photonic- 
crystal waveguides. Photonic-crystal waveguides are two- 
dimensional (2D) photonic-crystal slabs with a line de- 
fect wherein light is guided [8 j. They possess unique dis- 
persion relations, supporting slow-light propagation and 
enhanced light-matter interactions [9j [10] . 

With NSOM, one can measure the dispersion relation 
in these waveguides, map light pulses spatially and study 
slow-light propagation [3j [IT] [12] . It is also possible to 
measure the field patterns, which can be complicated be- 
cause of the multimodal nature of the structures. Spa- 
tial Fourier transforms are especially useful to analyze 
e.g., the dispersion [J or individual mode contributions 
[15] . For ballistic light propagation in photonic-crystal 
waveguides the detected field pattern is a superposition of 
Bloch modes determined by the symmetry of the waveg- 
uide. Ha et al. [14] recently proposed an algorithm that 
uses these symmetry conditions to extract Bloch modes 
from arbitrary measured field patterns [15] . In the optical 
domain the algorithm has so far been used to identify dis- 
persion relations in photonic-crystal waveguides [16] [17] . 
To date, however, 2D spatial patterns of Bloch modes 
at optical frequencies have not been obtained with this 
method. 

Here, we show the power of the Bloch mode recon- 
struction algorithm by extracting individual 2D mode 
patterns from phase-sensitive NSOM measurements on 
a GaAs photonic-crystal waveguide. We discuss in de- 
tail a specifically measured field pattern for which the 
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Figure 1: (color online) (a) Top view of a photonic-crystal 
waveguide that is divided in N unit cells (rectangles) for the 
reconstruction algorithm, (b) Calculated bandstructure for a 
GaAs photonic-crystal waveguide showing both TE-like (red) 
and TM-like (blue) guided modes. The gray area marks the 
light cone, LL is the light line, the blue (pink) area marks TM- 
(TE)-like slab modes, the purple area marks both TE- and 
TM-like slab modes. The labeled arrows mark the different 
modes considered here. 



Bloch modes are reconstructed. We apply this algorithm 
to study the behavior of the spatial width of the lowest 
frequency TE-like guided mode in the 2D band gap for 
TE-polarized light as a function of the wavevector. 



SAMPLES AND METHODS 

Figure []Ja) illustrates the top view of the waveguide 
studied here. It consists of a GaAs photonic-crystal slab 
with holes forming a triangular lattice with pitch size a = 
240 ± 10 nm, normalized hole radius ^ = 0.309 ± 0.002, 
waveguide length of approximately 1 mm and slab thick- 
ness h = 160 ± 10 nm. A single row of missing holes 
forms a Wl waveguide. Light is guided in the ^-direction. 
Each numbered rectangle represents a unit cell. Fig. 
[IJfr) shows the calculated bandstructure along x for the 
photonic-crystal waveguide surrounded by air using a 
plane wave expansion [18] (assuming a constant refrac- 
tive index of nc a As = 3.56 and thickness ^ = 0.67). The 



2 



black diagonal line represents the light line, correspond- 
ing to light propagation in air. Modes above the light 
line couple to modes outside the waveguide, are there- 
fore lossy [8 and are not considered here. The blue 
and pink areas mark a continuum of modes propagat- 
ing in the surrounding photonic crystal for TM- and TE- 
polarized light, respectively. The blue and red bands 
describe modes for, respectively, TM- and TE-polarized 
light that are guided by the line defect. Here we concen- 
trate on mode A that is TE-polarized, and modes B and 
C that are TM-polarized. 

A continuous- wave diode laser (Toptica DL pro 940) is 
used with a linewidth of 100 kHz and an emission wave- 
length between 907 — 990 nm, corresponding to a reduced 
frequency in the range 0.24 ... 0.26 a/ A. Light is side- 
coupled on a cleaved end-facet of the 1 mm long waveg- 
uide with a high-NA glass objective (NA=0.55). The 
incident light is linearly polarized with an angle of about 
45° with respect to the normal of the waveguide to ex- 
cite both TE- and TM-like modes. The field pattern is 
collected approximately 100 |Ltm away from the coupling 
facet using an aluminum coated fiber tip with an aper- 
ture of 160 ± 10 nm. We perform phase-sensitive NSOM 
using heterodyne detection. Detailed descriptions of a 
similar setup are presented in Ref. [TTJ 

The propagating modes considered here represent 
eigenmodes of the crystal and can therefore be repre- 
sented by Bloch modes. A 2D Bloch mode propagating 
in the ^-direction at position r = [x,y] and frequency uj 
is described by \I/ m (r,u;) = ^ m (r, cj) exp (z/c m ^) . Here 
%l) m (Y,ijj) is an envelope that is periodic with the lat- 
tice and satisfies ^ m (r,cj) = i/j m (r + ax,u), m labels 
the Bloch mode, and k m is the corresponding normalized 
Bloch vector (we consider normalized wavevectors only). 
We restrict ourselves to propagating modes (k G M). It is 
assumed that the measured field pattern <3>(r,co>) can be 
described by a superposition of M Bloch- modes \P m (r, uj) 
with amplitude a m and one overall residual e(r,uj): 
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Figure 2: (color online) (a) Measured amplitude for a 
photonic-crystal waveguide at 931.6 =b 0.1 nm. (b) Fitted 
amplitude using 7 reconstructed Bloch modes over 30 unit 
cells, (c) Residual amplitude, (d) Amplitude coefficients of 
the spatial Fourier transforms of the measured near-field pat- 
tern. Labels correspond to those in Fig. |TJ6) . 



suits in the field pattern A m = a m / m (r, uj) and the Bloch 
vector km for m = 1 . . . M. We compare the extracted k m 
with wavevectors determined from spatial Fourier trans- 
forms to confirm the accuracy of the algorithm. 



RESULTS 



M 

$(r, uj) = a m^m(r, uj) + e(r, uj) (1) 

rn=l 

The residual £(r, gj) describes measured field patterns 
that cannot be described by the M Bloch modes, such 
as non-guided modes, but also accounts for experimental 
artifacts and noise. 

In the reconstruction algorithm [HI [17] a section of the 
waveguide is separated into N unit cells, see Fig. [TJa). 
The algorithm uses the property that each ip m (r ,uj) is 
periodic in the photonic-crystal lattice, and requires us 
to analyze the measured field for each unit cell U n (r\uj), 
with r' the coordinate within one unit cell. The measured 
<£(r, uj) can be fitted with a series of Bloch modes using a 
least squares optimization that minimizes the functional 
W = J \s(r,uj)\ 2 dr/ j |<£(r, uj)\ 2 dr. This procedure re- 



Figure [2|a) shows the measured field amplitude 
|3>(r, tv)| for a waveguide section of 30 unit cells (7.2 ±0.5 
|Ltm) at a wavelength of A = 931.6 ± 0.1 nm. A beating 
pattern with a period of 5.7 ± 0.2 unit cells reveals that 
multiple modes are involved in the spatial pattern. To 
ensure that each U n (r\uj) describes precisely one unit 
cell, the original measurement was resampled on a differ- 
ent grid. Figure |5|6) shows the reconstructed amplitude 

M 

| ^m( r 9 ^)| with M = 7 Bloch modes, which is in 

m=l 

excellent agreement with the measured amplitude. We 
have chosen M = 7 because we expect from Fig. [Tj6) 
a forward propagating mode corresponding to the light 
line and for modes A, B and C both forward and back- 
ward propagation, hence M = 1 + (2 x 3) = 7. The low 
values for the functional W = 0.0133 demonstrates that 
indeed $(r, cv) is well described by the superposition of 7 
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Table I: Comparison between fitted wavevectors from Bloch- 
mode reconstruction (km), and obtained wavevectors from 
spatial Fourier transforms (&sft). The first column labels the 
Bloch modes. The second column gives the fitted k m . The 
third column gives a measure how strongly present a mode 
is. The fourth column gives the /csft, where the superscript 
F marks the fundamental wavevector. In the fifth column we 
identify the modes from the calculated bandstructure in Fig. 
|Tj6) , the propagation direction, where +( — ) corresponds to 
the positive (negative) x-direction, and polarization. 



k m ( 2 ™) 



feFT (- 



Label 



0.251(1) 0.0073 

0.330(3) 0.0117 

0.379(2) 0.0086 

0.488(2) 0.6714 

-0.331(2) 0.2202 

-0.367(6) 0.0014 

-0.486(2) 0.0681 



0.251(3) 
-0.668(3) F , 0.331(3) 

0.370(3) 
-0.516(3), 0.489(3) F 
-0.332(3), 0.668(3) F 

-0.370(3) 
-0.489(3) F , 0.516(3) 



+LL 
-A, TE 
+B, TM 
+C, TM 
+A, TE 
-B, TM 
-C, TM 



Bloch modes. This conclusion is confirmed by the abso- 
lute residuals \s(r,uj) \ plotted in Fig. [2jc) . The fitted k m 
are presented in the second column of Tab. |TJ The third 
column describes the relative contribution of each mode 
as c m = I f A m (r,uj)$(r,uj)dr\/ J |$(r, uj)\ 2 dr. Note that 
the 7 contributions plus that of the residual add up to 
unity. The fifth column describes which modes of Fig. 
[Tj6) correspond to fc m , the propagation direction and po- 
larization. We observe mainly the forward propagating 
TE-like mode A (m = 5) and the forward propagating 
TM-like mode C (m = 4). The errors in k m are esti- 
mated by varying the grid element size and allowing for 
a relative increase of AW by maximum 10%; within this 
range the mode patterns A^t^uj) and A 5 (r,c<j) do not 
change noticeably. 

Next, the fitted k m are compared with wavevectors 
determined from the spatial Fourier transforms shown in 
Fig. (2jcf) (&sft)- A Fourier transform in the x -direction 
was made for each line parallel to the waveguide over a 
range of 35.4 ± 0.9 |im, which includes the range shown 
in Fig. [2|a). For comparison k m and /csft are listed 
in Tab. [ij showing an excellent agreement. The spa- 
tial Fourier transforms show for modes A and C higher 
Bloch harmonics. Both the fundamental /csft and the 
observed higher Bloch harmonics are listed in the table. 
In Fig. ^d) the modes from Fig. [TJfr) are identified. The 
amplitude coefficients confirm that we detect mainly the 
forward propagating TE-like mode A (&sft = —0.332 ± 
0.003, 0.668 ± 0.003) and the forward propagating TM- 
like mode C (fc S FT = -0.516 ± 0.003, 0.489 ± 0.003). 

We have demonstrated that the forward propagating 
TE-like mode A and the forward propagating TM-like 
mode C are the most prominent Bloch modes present 
in the data of Fig. |5J Figure ^a) shows the amplitude 
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Figure 3: (color online) (a) Reconstructed amplitude for Fig. 
|5Ja) using Bloch modes m = 4 and m = 5 only, (b — c) Am- 
plitude, imaginary part and real part for both Bloch modes. 
(d — e) Comparison between calculated \E\ (left) and recon- 
structed \A\ (right). 



when only these two modes are taken into account for the 
reconstruction. A very good agreement is observed with 
&(r,u) of Fig. [5|a). Especially the diagonal beats are 
well reproduced. The difference wavevector of the two 
modes corresponds to a beating period of (5.5 ± 0.2)a. 
The beating pattern of two orthogonal modes is the re- 
sult of quasi-interference; the NSOM tip thereby projects 
both orthogonal polarizations on a detection basis where 
these modes interfere [19 . Figure [3^6) shows the ampli- 
tude, the real part and the imaginary part of the recon- 
structed TE-like Bloch mode A with k 5 = -0.331 ±0.002. 
The mode profile is symmetric in the ^-direction about 
the center of the waveguide. Figure ^d) shows the cal- 
culated [18] time-averaged amplitude (\E\) (left) and the 
measured amplitude \A\ for approximately 3 unit cells. 
Both show a similar pattern. Figure [3^c) shows the am- 
plitude, the real part and the imaginary part of the recon- 
structed TM-like Bloch mode C with /c 4 = 0.487 ±0.003. 
Figure [3^e) shows the calculated (\E\) (left) and the mea- 
sured amplitude \A\ for approximately 3 unit cells (right). 
For mode C the agreement is poor, likely because the 
near-field tip has a low response to E z and a non-trivial 
response to E x and E y: see Ref. 13, in addition to its 
finite resolution. The reconstructed field patterns are an 
approximation of the Bloch modes propagating in the 
system and moreover, they are not necessarily orthogo- 
nal because of quasi-interference. Although the effect of 
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Figure 4: (color online) (a) Measured normalized transverse 
mode profile for k x = 0.285 (black), k x = 0.332 (red dashed), 
and k x = 0.375 (blue), (b) Calculated normalized transverse 
mode profile for k x = 0.280 (black), k x = 0.330 (red dashed), 
and k x = 0.380 (blue), (c) Determined (black) and calculated 
width w (red) versus longitudinal wavevector k x for the TE- 
like guided mode. Inset: measured (symbols) and calculated 
dispersion (red). 



the tip is far from straightforward ([5, 20, 21 J, we antic- 
ipate that the comparison between calculated modes of 
an optical system and reconstructed modes could lead to 
methods to deduce the response function of a near-field 
tip. 

Next, we demonstrate the power of the reconstruction 
algorithm by studying the transverse behavior of TE-like 
mode (A) versus wavevector. We consider the forward 
propagating TE-like mode (fc 5 = -0.331, 0.669 in Fig. [3| 
and take its reduced wavevector (0 < k x < 0.5) to com- 
pare directly with the folded bandstructure of Fig. [IJfr). 
We have measured field patterns between A = 907 — 944 
nm, and apply the reconstruction algorithm to determine 
A(r, uj) and k x for this mode at each A. We have selected 
3>(r, a;) where this TE-like mode is prominently present 
in spatial Fourier transforms. The inset in Fig. [4jc) 
shows the fitted k x versus reduced frequency j. 

In order to concentrate on the transverse behavior, we 

a 

define the transverse mode profile J |A(r, uS)\dx. Fig- 

o 

ure[4ja) shows the measured normalized transverse mode 
profile for 3 different wavevectors. At k x = 0.285 (black) 



a transverse mode profile is apparent that can be mainly 
described by one prominent maximum at ^ = that is 
slightly asymmetric, describing light guided in the line 
defect. Additional side lobes are observed at - = —0.9 

a 

and at ^ = 1.4, representing light extending into the sur- 
rounding photonic crystal. At k x = 0.332 (red dashed) 
we observe a central maximum at - =0, and the con- 
tributions of the side lobes become bigger. Also note 
the new peaks observed at ^ = —3.2,-2.2,-1.5, and 
2.8. When the wavevector is increased, the relative 
contributions of these additional peaks increase. At 
k x = 0.375 the central peak is still present, and the 
surrounding peaks have grown. Qualitatively, the mea- 
sured transverse mode profiles correspond with the calcu- 
lated ones shown in Fig. |4^6), which were obtained from 
the time-averaged amplitude of the total electric field 

a 

J \E(r, uS)\dx. The maxima and minima occur at approx- 
o 

imately the same locations and the width w of the central 
maximum is growing with increasing k x . Not all features 
are resolved of the calculated transverse mode profile in 
our measurements. For example, the measured relative 
amplitude of the central maximum compared with the 
additional maxima differs from the calculations. 

The central maxima of the transverse mode profiles are 
fitted with a Gaussian with width w. Figure [4jc) shows 
the reduced width ^ versus wavevector (black symbols, 
bars represent 95% convergence intervals). The red sym- 
bols interpolated by the dashed line represent ^ deter- 
mined from the calculated transverse mode profile. The 
measured w increases with k x . For k x < 0.34 the mea- 
sured w matches the calculated w well. For k x > 0.34 
the measured w becomes larger than the calculated w. 
We attribute this to the finite resolution of the near-field 
tip and to its response function. For most considered 
/^-values, however, the measured w/a matches the cal- 
culated with a 10 % accuracy. 



SUMMARY 

In conclusion, we have implemented an algorithm pro- 
posed by Ha et al. to extract Bloch modes from near-field 
measurements on a photonic-crystal waveguide. The ex- 
tracted wavevectors are in very good agreement with the 
wavevectors determined from spatial Fourier transforms. 
We have studied two extracted Bloch modes to explain 
the observed near-field pattern. We also have studied 
how the width of a selected mode changes with wavevec- 
tor and find good agreement with calculations. We antic- 
ipate that this algorithm can be used to filter states that 
cannot be described by propagating Bloch modes, such 
as Anderson-localized states observed in the slow-light 
regime [22H26] . 
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